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Abstract 

G ■ Maximally random jammed states of hard spheres are prototypical glasses. We study the small 



wavenumber k behavior of the structure factor S{k) of overcompressed million-sphere packings as 



CN ! a function of density up to the jammed state. We find both a precursor to the glassy jammed 

> 
0^ ' state evident long before the jamming density is reached and an associated growing length scale, 



extracted from the volume integral of the direct correlation function c{r)^ which diverges at the 



00 

en . "critical" jammed state. We also define a nonequilibrium index X and use it to demonstrate that 

CN I the packings studied are intrinsically nonequilibrium in nature well before the critical point is 



reached. 



A sufficiently rapid quench of a liquid from above its freezing temperature into a super- 
cooled regime can avoid crystal nucleation to produce a glass with a relaxation time that is 
much larger than experimental time scales, resulting in an amorphous state (without long- 
range order) that is simultaneously rigid jl| . The underlying physics of the glass transition 
is one of the most fascinating open questions in materials science and condensed-matter 
physics. Many conundrums remain, including whether the growing relaxation times under 
supercooling have accompanying growing structural length scales. Essentially, two opposing 
explanations have emerged to address this question. One asserts that a static structural 
length scale does not exist and identifies growing dynamical length scales |2rl4[. The other 
contends that there is a static growing length scale of thermodynamic origin 5|, l6|. In this 
paper, we present both theoretical and computational results that support an alternative 
view, namely, the existence of a growing static length scale but one that is intrinsically 
nonequilibrium in nature. 

Our model systems are disordered packings of identical spheres with densities between 
the freezing transition and the so-called maximally random jammed (MRJ) state 7|. The 
MR J state under the strict-jamming constraint is a prototypical glass in that it lacks any 
long-range order but is perfectly rigid (the elastic moduli are indeed unbounded) |8|, |9|. This 
endows such packings with special attributes. For example, MRJ packings are hyperuniform 



11| {i.e., infinite wavelength density fluctuations vanish) with a structure factor S{k) 



that tends to zero linearly in the wavenumber k, implying quasi-long-ranged negative pair 
correlations (anticorrelations) decaying as a power law [UJ]. This large-scale property is 
markedly different from typical liquids in equilibrium, which tend to exhibit more rapidly 
decaying pair correlations (including exponential decays). 

It has been theoretically shown that hyperuniform point distributions are at an "inverted" 
critical point, i.e., in contrast to normal fluid critical points, the direct correlation function 
c{r), rather than the total correlation function h{r) = (72 (r) — 1, with (72 (r) the pair corre- 
lation function, becomes long-ranged (d ecay ing more slowly than 1/r'^, where d is the space 



dimension and r is the radial distance) 



10| . This suggests that the volume integral of c{r) 



for overcompressed hard-sphere configurations 12| that follow Newtonian dynamics should 
grow as the MRJ state is approached and ultimately diverge at this inverted critical point. 
In this Letter, we show that this not only is the case but that the nonequilibrium signature of 
quasi-long-range anticorrelations, which we quantify via a nonequilibrium index X, emerges 



well before the jammed state is reached. This implies that the direct correlation function of 
a glass formed by supercooling a liquid (in which the molecules possess both repulsive and 
attractive interactions) provides a static growing length as its glass (jamming) transition is 
approached. Hence, our findings, based on nonequilibrium hard-sphere systems, could be 
applied more broadly to general glass formers, as we will elaborate below. 

The present study focuses on prototypical glassy states represented by three-dimensional 
identical nonoverlapping spheres. For this system, it has been shown that MRJ packings are 



characterized by quasi-long-ranged pair anticorrelations in which h{r) decays as — l/r*^ 11 |. 
This quasi-long-ranged behavior in r is equivalent to linear behavior in the structure factor 
S{k) near the origin with S{k) nonanalytic at A; = 0. The structure factor S{k) is defined 
in terms of the Fourier transform of h{r), S{k) = 1 + ph{k), with p the number density, 
equal to {Q/tt)(J) for unit-diameter hard spheres with packing fraction the fraction of space 
covered. Quasi-long-ranged anticorrelations have been shown to be present in the ground 
states of liquid helium 13| and noninteracting spin-polarized fermions [IJ], in the Harrison- 



Zeldovich spectrum of the early Universe 15|], and in MRJ states of Platonic solids |l6 |. 



Experimentally, the most accurate measurements of the structure factor of bulk amorphous 

l-k behavior of S{k) toward 5'(0) = that 
17j . Linear behavior near the origin in S{k) 



silicon have revealed a linear trend in the sma 

appears to be consistent with hyperuniformity 

for a hyperuniform system indicates an inverted critical point, characterized by long-range 

anticorrelations in a direct correlation function c{r) that decays as — 1/r^, where we define 

c{r) in terms of the Fourier transform of c{k), 

_m_ sjk) - 1 

'^^^ " s{k) - ps{k) ■ ^^^ 

More information about this inverted critical point, including critical exponents, can be 
found in Ref. |10 |. 



In this work, we study million-particle packings of identical spheres overcompressed at 
various rates fast enough to avoid the formation of crystallites. We employ an event driven 
friction-free Lubachevsky-Stillinger molecular dynamics protocol under periodic boundary 



conditions [18| in a cube of side length unity where initial sphere velocities are Maxwell- 
Boltzmann distributed such that mean energy per sphere (of mass unity) is fc^T = 1/2, with 
ks Boltzmann's constant and T temperature. Simulations are performed to overcompress 
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FIG. 1: (color online). Bottom curve: angularly-averaged S{k) calculated by DFT for a million- 
sphere packing with (j) = 0.9999999(/)c .... Note the linear behavior in S{k) as k ^>- 0. Top curve: 
average of 20 angularly- averaged S{k) calculated by DFT for independently generated million- 
sphere packings in the equilibrium liquid state with <j) = 0.45. Vertical bars represent an estimate 
of the standard deviation of the distribution of values for S{k) at each k over the 20 packings. 

spheres from various initial states (including equilibrium liquid) at densities below freezing 
(pf = 0.494 up to various percentages of the jamming fraction (pc- For each packing, (pc is 
dependent upon initial conditions and compression rate, but it varies no more than 5 x 10~^ 
for packings compressed at the same rate due to the very large system size. We study rates 
in this work corresponding to sphere diameter growth per unit time from F = 0.0007 to 
F = 0.03 181], where F = 0.0007 appears to be the slowest rate at which crystallites do not 
form. 

We closely examine the small wavenumber k behavior of the angularly-averaged structure 
factor S{k), where k is the magnitude of a wavevector k. We calculate S{k) by direct Fourier 
transform (DFT), omitting forward scattering, for A^ = 1,000,000 unit-diameter spheres in 
a periodic cube of length L replicated over all space, using S'(k) = (l/iV)| J2j=i e^^'^^\^, and 
angularly averaging over all k of equal magnitude k. The smallest nonzero wavenumber k 



calculated by this method is equal to ^Ti /L, which for 1, 000, 000 unit-diameter spheres at 
jamming is about 0.067. Figure [1] compares the behavior of S{k) for a) a single packing 
overcompressed from the equilibrium liquid at = 0.690c to = 0.99999990c and b) 
an average of 20 independently generated packings in the equilibrium liquid state with 
= 0.45 ~ 0.700c. 
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FIG. 2: (color online). Linear fits with at various percentages of jamming fraction 0c to the 
average of 20 angularly-averaged structure factors S{k) calculated by DFT from million-sphere 
packings compressed at a rate F = 0.03 from the equilibrium liquid. All fits achieve R-squared 
values of greater than 0.995. The curves in the figure from top to bottom correspond to the labels 
in the legend from top to bottom. Inset: average data and estimate of the standard deviation (over 
the 20 packings) for angularly- averaged S{k) for packings compressed to = 0.970c. 



For all compression rates and initial conditions studied, we identify a precursor to quasi- 
long-ranged pair anticorrelations manifest in linear behavior in S'(fc), i.e., S{k) = ak + b. 
For greater than about 0.920c, this behavior extends from the smallest k studied to about 
k ~ 0.4, whereas for 0/ < < 0.930c, a linear trend is clearly evident but does not extend 
to the smallest values of k for all compression rates studied. This behavior is illustrated in 
Fig. |2]for a wide range of densities |l9 |. 



Though plots of S{k) from individual packings show variability at small k around the 
linear trend attributable to the dynamics of particle movements as the spheres are com- 
pressed, the average of several packings always illustrates a distinctly linear trend. This 
trend additionally appears for packings compressed at various rates from a random sequen- 
tial addition (RSA) state |2Q|, and the trend appears for packings compressed at fast but 
changing rates. Taken together, these observations indicate that the linear trend is easily 
reproducible across different compression rates and initial conditions, even nonequilibrium 
conditions like RSA, so long as the compression rate is fast enough such that crystallites do 
not form. 

As — )■ 0c! the the values of c{k) (Eq. ([T])) near k = begin to diverge. In particular, in 
the limit as — > 0c for a hyperuniform packing, c(0) — )■ — oo. This divergent behavior can 
clearly be seen in Fig. [3l which are plots of c{k) calculated using Eq. ([1]) and the fits in Fig. 

m 

The length scale ^dcf = (^c(O))^/^ grows continuously with packing fraction and diverges 
as the hyperuniform state is approached, indicating a long-ranged direct correlation function 
c{r) that decays asymptotically proportional to — 1/r^. The growth in C,dcf signals the 
incipient rigidity associated with the formation of nearly jammed sphere-contact networks, 
and C,DCF diverges when the packing becomes completely rigid at the inverted critical point 
= 0c. The presence of this inverted critical point suggests that renormalization group 
methods might be fruitfully applied to study the behavior of c{r) in this and similar glassy 
systems. 



In the mean spherical approximations 2]J], c{r) can be thought of as the negative of an 



effective pair potential Veff{r), that is to say, —/3veff{r) = c{r) with (3 = l/fc^T, which 
is asymptotically exact as r — )■ oo. Following this interpretation, MRJ packings exhibit 
long-ranged repulsions fe//(r) that are asymptotically proportional to 1/r^, which can be 
thought of as a generalized Coulombic interaction 1/r" with n = 2 instead of n = 1, that 
drive the system to have no infinite wavelength density fluctuations, i.e., to be hyperuniform. 
In comparison, one-component plasmas in equilibrium interacting with a Coulomb potential 
(e.g., v{r) ~ 1/r in 3D) are hyperuniform id. Il4j. |22|. 

Other diverging length scales can also be identified, e.g., the inverse of the first point 
k in S{k) where linear behavior becomes dominant, which we term ^sf- This length scale 
is visually evident in S{k) for the systems studied beginning at just above 0j, and it 
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FIG. 3: (color online). Fits of angularly-averaged c{k) calculated from the fits in Fig. [2j Note 
that c(0) appears to diverge to — oo as — )• cpc- Inset: average data and estimate of the standard 
deviation (over 20 packings) for c{k) for packings overcompressed to (/> = 0.97(/)c. For this data set, 
geometric {fig = (ni=i ^j)^ ) with fig the mean and Xi the data points) and arithmetic averages 
were nearly identical. 



can be shown 23|] to correspond roughly to the point in sphere packings where — 1/r^ pair 
anticorrelations in h{r) cease to be dominant and are supplanted by faster decay. This 
length scale C,sf grows far more quickly than ^dcf] in particular, it is already greater than 
the size L of the million-sphere systems studied once = 0.930c- 

The identification of static growing length scales that can be extracted from pair infor- 
mation is a novel and unexpected finding. Though diverging time scales have long been 
evident in glasses, quickly growing length scales have been difficult to find. Recent studies 
have highlighted the presence of growing length scales using both static and dynamic four- 
point correlation functions; however, none of these length scales have been shown to grow 



very large, and it has been suggested that standard pair correlation functions may 
information necessary to describe the static growing length scales present in glasses 



ack the 



24 



25|. 



Our results demonstrate otherwise. By refocusing investigations to examining the small-fc 



behaviors of c{k) and S{k), which are accessible to us via very large million- sphere packings, 
we are able to determine the large-r behavior of c{r) and h{r) (which is difficult to do in 
real space) to extract growing length scales. 

Though there are diverging length and time scales evident in the structure factors of these 
prototypical glasses, the packings are far from equilibrium. In particular, the compressibility 
equation relating the isothermal compressibility kt = {l/(f)){dp/d(j))T, with p the pressure, 
to infinite wavelength density fiuctuations, 

5(0) = pknTKT, (2) 

which holds for systems in thermal equilibrium, does not hold for the systems studied. To 
investigate this concept further, we make use of an expression for the reduced pressure 
p/pksT that is valid continuously from the freezing point (pf all the way up to the jamming 
density (pc as obtained from the nearest-neighbor conditional pair distribution function G{r) 
|20|. Specifically, 

G(oo) = p/pksT = 1 + 4</>^^(l)^^, (3) 

where 5^/(1) is the value of (72(1) for an equilibrium liquid at 0/. From Eq. (J3]), we can derive 
the right hand side of Eq. ([2]), 

pkeTKT = -j:^ T-TT- —T- (4) 

cC// 1 0C-0/ 



Relation ([3]) has been shown to fit closely to data 26|], and our results, displayed in Fig. HJ 
support this finding. Values of pkBTux calculated from Eq. (jlj) also fit the data. In Fig. HJ 
we calculate kt by rescaling velocities after each particle has undergone two collisions so that 
temperature remains constant. The expression pksTKT is computed once every time parti- 
cles have undergone 20 collisions each from the relation pksTnx = {l/(j)){A(j)/A{p/pkBT)), 
with p the pressure, and the average over ten calculations is taken. The average of these 
values over five packings is reported as the measured isothermal compressibility. Though 
reduced pressure, and therefore kt, should change if spheres are allowed to relax (collide 
with no compression), in our simulations this change was undetectable within the error even 
for values of pksTuT after relaxation times of 100, 000 collisions per sphere at = 0.930c- 
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FIG. 4: (color online). Average over 20 packings compressed at rate T = 0.03 of reduced pressure 
at various percentages of (j)c, compared to the calculated values of reduced pressure from Eq. ([3]). 
Inset: average values of pksTKT at each value of ^ for four packings, compared to the calculated 
values using Eq. @. The standard deviation of the measurement is about the size of the circle 
displayed for (j) = 0.93</»c and much smaller than the circles for the other three measurements. 



With Eq. (jlj) in mind, we define 



X 



S{0) 
pksTKT 



(5) 



as a nonequilibrium index to quantify the degree to which the systems under study devi- 
ate from thermal equihbrium X = [cf. Eq. (l2])]. Figure |5] shows a plot of X vs. 0, 
demonstrating that X diverges as — )■ (pc- 

Using Eq. (jl]) and fitting a hnear trend to values of S'(O) with from Fig. [2] (inset of Fig. 
|5]), the behavior of X as — )■ 0c for these systems can be calculated. Though both pksTKx 
and 5(0) approach zero as — )■ 0c, X still diverges with a pole of order one at = 0c. This 
is an interesting result, and it strongly indicates that the jammed glassy state for this model 
is fundamentally nonequilibrium in nature, though at this time it is unclear precisely how 



X will behave near 0c for other glassy states and hyperuniform systems. 
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FIG. 5: (color online). Plot of X with phsTuT calculated from Eq. (j4]) and 5(0) taken from the 
linear fits in Fig. [21 compared to where 5(0) is calculated from a linear fit to the intercepts 5(0) 
in Fig. [2j Inset: linear fit to the intercepts 5(0) in Fig. [2l 

We have shown that for a variety of mono disperse hard-sphere systems compressed at 
various rates from various initial configurations that there is a precursor to the jammed 
glassy state evident far below the jamming packing fraction 0c- This precursor appears as 
linear behavior in the structure factor S{k) near A; = associated with a quickly growing 
lengthscale ^sf-, and for hyperuniform systems is indicative of the onset of an inverted critical 
point at = 0c associated with a diverging length scale ^dcf and long-ranged power law 
decay of the direct correlation function c{r). Our identification of these length scales is a 
unique finding that demonstrates the onset of the jammed glassy state is detectable using 
standard pair information [271 ]. 

Due to the early onset and robustness of the jamming precursor and the dynamics em- 
ployed in our model of a prototypical glass, we expect that our results are broadly applicable 
to glass forming molecular systems {e.g., metallic glasses, network glasses, etc.). Hard-sphere 



10 



packings nearjamniing are known to provide excellent structural models of glassy molecular 



states 



20 



Endowed with Newtonian dynamics and quickly compressed, hard-sphere 
models are driven by strong pair repulsion and free-volume dynamics, which are the salient 
drivers in their molecular counterparts. These drivers lead to diverging elastic moduli and 
timescales just as sudden macroscopic rigidity and diverging timescales appear in molecular 
glass formers. Therefore, the precursor to the jammed glassy state and the diverging length 
scales that we have shown to be present in hard-sphere glasses are applicable broadly to 
other glass formers [291]; prehminary results indicate that ^dcf and ^sf grow as a function 
of inverse temperature in a supercooled model molecular glass former. 

Despite the presence of diverging length scales, we have shown via a nonequilibrium index 
X that hard-sphere glasses are intrinsically nonequilibrium in nature, where X appears to 
diverge as — )■ (pc- Though it is currently unclear if X will diverge similarly in supercooled 
model molecular glass formers, we expect that it will generally increase with decreases in 
temperature as glassy states are approached. In future work, we will also investigate the 
behavior of X as a function of temperature in supercooled model molecular glasses. 

This work was supported by the NSF under award number DMR-0820341. 
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